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Abstract 

Quenched chiral perturbation theory is used to compute the first finite volume correction to the 
chiral condensate. The correction diverges logarithmically with the four-volume V. We point out 
that with dynamical quarks one can obtain both the chiral condensate and the pion decay constant 
from the distributions of the lowest Dirac operator eigenvalues. 
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1 Introduction 



Finite-size effects in lattice gauge theory are usually considered undesireable features of numerical 
simulations. There are exceptions: near critical points, for instance, finite-size scaling applies the 
Renormalization Group in finite volumes to extract critical indices of the infinite volume system. 
Also away from critical points one can use knowledge of the finite volume behavior to learn about 
the infinite volume theory. When this is the case, one is turning the apparent disadvantage of being 
restricted to small volumes into a distinct advantage. 

Ordinarily finite-volume computations in gauge theories still involve the full complexity of such 
strongly coupled theories. A program to overcome this obstacle was initiated several years ago by 
Gasser and Leutwyler |^. Let the finite (euclidean) 4-volume be of size V = LxLyLzLt{Lt = l/T). 
The idea is that for Lj » 1/Aqcd the euclidean partition function (or appropriate generating func- 
tional) is dominated by the lowest-mass excitations. If chiral symmetry is spontaneously broken in 
the infinite volume limit these light particles are the pseudo-Goldstone bosons, whose interactions 
vanish at zero momentum transfer. The theory governing these excitations is a chiral Lagrangian. 
Although non-renormalizable, this theory can be treated perturbatively in the sense of an effective 
Lagrangian. Combined with the description of the fundamental degrees of freedom in the underlying 
gauge theories, this allows for precise finite- volume computations in these otherwise strongly coupled 
theories. 

There are two finite- volume regimes to consider. The first is the one connected with the usual infinite 
volume limit of the theory. If M denotes the mass of the pseudo-Goldstone boson, one requires 1 /M <C 
Li so that the particle "fits well inside the box", and the finite volume y is a minor perturbation. 
In this regime usual chiral perturbation theory is applicable. The other regime is 1/M ^ L^: if the 
volume is much larger than the scale set by the gauge theory, Lj ^> 1/Aqcd, the euclidean partition 
function is still dominated by the lightest excitations, but ordinary chiral perturbation theory fails. In 
this region a different expansion must be applied |^]. This paper concerns (partially) quenched theory 
in the latter regime. We shall also comment on how a known result from the unquenched theory can 
be used to extract the pion decay constant F in an unusual, but quite simple way. 

It is worthwhile stressing that although the finite volumes which are being probed here are unphys- 
ical in the sense that the pion correlation length exceeds the boundaries of the finite volume, the 
observables one can extract in this way are infinite-volume quantities. It is precisely the purpose to 
use unphysical finite-volume expectation values to extract these infinite-volume quantities directly, 
without any numerical extrapolation to the infinite- volume limit. 

2 The Chiral Condensate 

We begin with a brief review of the Feynman rules for partially quenched chiral perturbation theory 
based on the replica method Denoting the number of ordinary physical ("sea") quarks by Nf, one 
adds (in its simplest form) "valence" quarks of common mass niy to theory, and takes the limit 
A'^ — > in the end. The starting point is then an effective theory based on the coset of 

\J{Nf + N,) X \J{Nf + N,) ^ \JiNf + Ny) , 

including the flavor singlet field. The Lagrangian can be parametrized 

C = ^Tr{d^U{x)d'^U\x)) - ^TrM{U{x) + U\x)) + + ^d,M^)d^Mx) • (1) 
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when keeping only the usual leading terms in the small momentum expansion. The field U(x) = 
ex.p[i^/2^(x)/F] G U(A^/ + N^) lives on the Goldstone manifold, and we have separated out the flavor 
singlet combination ^q{x) = Tr<I>(x). For this reason it is convenient to work in the quark basis 0] 
where ~ 'ipiipj- Let M^j = {rrii + mj)'S / F'^ , and 

Hp') = i + ^^^^li^^ + i:-r^l • (2) 

The momentum-space propagator for the off-diagonal mesons ~ ipiipj, i ^ j is then read off to be 

while the corresponding propagator for the diagonal combination $jj ~ ipii^i is 

^ (2. ^ Sjj (mg + ap^)/Nc 

With these rules it is straightforward to take the partially quenched limit of chiral perturbation theory: 
one simply sends A'^j, ^ at the end of the calculation. It is simple to show that it is equivalent 
to partially quenched chiral perturbation theory based on the supersymmetric formulation @, ^, ^, 
but we find the replica method easier to use. It is straightforward to extend the method to a more 
elaborate chiral Lagrangian involving higher order terms in the momentum expansion. Here we restrict 
ourselves to the 0{p'^) Lagrangian, supplemented with the lowest order <I>o-terms (|l|). 

The conventional regime in which to apply (partially quenched) chiral perturbation is the one where the 
four-volume V well encloses the pseudo-Goldstone bosons of masses Mij, i.e., L » l/Mij. Although 
the four-volume V then can be kept finite, the correlation length does not reach from boundary 
to boundary, and the theory looks like a mildly perturbed infinite-volume theory. In this regime 
chiral perturbation theory can be arranged as a systematic expansion, roughly speaking in the small 
momentum p. This is clearly an obvious regime of interest for lattice gauge theory. 

Actually, as pointed out first by Gasser and Leutwyler [|^, another finite- volume regime can also be 
used to extract physical observables. This is the extreme finite- volume regime where L <C l/Mij. 
Although the pseudo-Goldstone bosons here apparently violate the bound which one would normally 
impose on the finite-size system, one can still perform analytical calculations here, and thus extract 
physical quantities by comparison with Monte Garlo data. This regime is governed by a perturbative 
theory of the chiral Lagrangian which is sometimes called the e-expansion (not to be confused with 
an expansion away from critical dimensions). Let us first introduce the mean length scale 

L = . (5) 

The counting of orders of e is defined by fixing 1/L = 0{e). The non-zero momentum modes thus 
count as p ~ 0{e) too, while rrii ~ ©(e^) (and hence Mij ~ 0{€'^); fii = uiiTiV is of order 1). A 
systematic expansion can be performed in which one is working to a given fixed order in e Q. 

In the above regime the leading term of the chiral Lagrangian is given by the zero momentum mode 
of the field U{x). The method of ref. H] treats this case by means of a collective field technique. To 
this aim, introduce a new variable ^(x) containing non-zero momentum modes only, and write 

U{x) = ue*^«(")/^n , (6) 
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where u € \J{Nf + N^) is a constant collective field corresponding to the zero momentum modes of 



the pseudo-Goldstone bosons. Let Uq '■ 



and introduce the flavor singlet 



with Ui 







exp 



[iV2^o/F] . 



(7) 



and E{x) = Tr^(x). The change of variables (|6|) is next inserted into the path integral. The Jacobian is 
non-trivial, but to lowest order it contains a term that only shift the vacuum energy, without affecting 
the effective Lagrangian calculation below. We then find that the action to lowest order becomes 



dx C 



dx Tr 
^Tr 



-dMx)d''c{x) + ^ 



M{Uo + U^o)(V 



F2 



- E 

dx ^{xf 



a 
2N, 



■df,E{x)d^'Eix) 



(8) 



It is convenient to add to this action a mass term for the field (,{x), and only taking this mass to 
zero at the end of the calculations. Because of the Gaussian damping term for the non-zero modes, 
^(x) counts as being of order e. For the singlet field E(x) the counting is not as straightforward. 
The kinetic energy term for E(x) suppresses fluctuations in E{x) beyond 0(e), so that the mass 
term for E{x) is of order e(mo)^e~^. The parameter mo is here a priori free as far as counting in 
e is concerned, and we thus have exactly the same counting ambiguity as in conventional partially 
quenched chiral perturbation theory [Q, ||: to any given order in the momentum expansion we should 
in principle include all orders of at each order in M^j. What we need to know here is that we 
can take tuq/ [NcF"^) to be parametrically small, and thus perform a simultaneous expansion in e and 
ml/{NcF'^). This is what we will do below. 

Integrating out the fluctuation field ^(x) to one-loop order, we obtain the one-loop improved effective 
Lagrangian for the zero momentum modes. We are particularly interested in the first correction 
to the chiral condensate S. What enters repeatedly in what follows is the pseudo-Goldstone boson 
propagator, evaluated at zero, 

,2^ - ■ X- ■ (9) 



A(Mfi)-fi:^. 



and the corresponding quantity for the modes without zero momentum, 



r2^ ^ 1 



- V- 



(10) 



The two are hence related by A(M2.) = (M^ F)"^ + A{M^^). 



It is instructive to first rederive the known finite-volume correction for the unquenched theory 
the present formulation. We thus evaluate the one- loop ^-integral saturation of 



m 



2F2 



Tr 



where for convenience we have taken the Nf dynamical fermions to be of equal mass (and the 
pseudo-Goldstone bosons thus have common masses Mss = 2ms S/F^). Taking the = Q expressions 
from above, we get^ 



dx 5C 



^Tr(C/o + Ul) {Nf - l)K{Ml) + - ^ 'G{MI 



^The coupling between Ho and 'S.{x) is treated perturbatively. It contributes only to the next order. 
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Tr(C/o + Ul) {Nj - l)K{Ml) + h{Ml) - -J2 



[11) 



This differs from the expression quoted in ref. lH because we have here included the effect of the flavor 
singlet field. Due to the anomaly, the mass of this field stays finite in the full theory as the quark 
masses are decreased, while the pseudo-Goldstone bosons approach zero mass. This corresponds to 
taking ttt-q — > oo in the above expression, which after setting the argument of the propagator to zero 
gives 



dx 5C 



2F2 



TV(f/o + [/, 



N 



f 



1 



N 



f 



-A(0) 



(12) 



and we hence recover the result of ref. H. 



We now turn to the fully quenched case. This simply corresponds to setting Nf = in the Feynman 
rules given above. We evaluate 



lim 



1 I m^iS 



2F2 



Tr 



(f/o + Vl) J dx 



which is a finite quantity. In complete analogy with eq. (Ill]), this results in 



lim 



1 Vm^T; 
JV^-^o iV^ 2F2 

1 VnivTi 
hm —r- 

N^^o 2F2 

1 Vm„Ti 



1 



Tr(C/o + C/J) <j (iV. - l)A{Ml) + - ^ 'G(m2J 



Tr(^o + t/J) (N, - 1)A(M2J + A{MI) - -J2 



— lim 



N^^oN^ 2F2 



Tr(C/o + C/t)-^ 



(13) 



after doing the integration over the non-zero modes only. As expected, this involves the quenched 
double pole |, |]. We can rewrite the correction as 



1 Vm^T, 



lim 



TV(C/o + U^o) {«A(m2J - {ml - oM^JSm^XM^j} 



Let us summarize what this correction implies. To zero'th order, the mass-dependent quenched chiral 
condensate in finite volume V is simply given by a group integral. We introduce the partition function 
at zero 0- angle 



dUo exp 



2Nr 



(14) 



It follows from the above analysis that this is the zero'th order approximation of the effective euclidean 
partition function of QCD with A^^, mass-degenerate (light) quarks in the finite-volume regime V ^ 
I/Aqqjj and <^ l/V |^], including the effects of a possibly light flavor singlet. Use has been 
made of the fact that the first term in the effective Lagrangian depends only on the combination 
fiy = m^VE. Formally the quenched chiral condensate then follows from 



(15) 



While it may be difficult to evaluate this expression explicitly, both low-mass and large-mass series 
expansions have been obtained in sectors of fixed topological gauge field charge & P]. 
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To one-loop order in the fluctuations of the non-zero momentum modes we see that this result is 
modified in a simple way. One can use the same effective partition function ( [T^ ) after replacing fiy by 



i 1 + 



aA(M2j - {ml - aM^ja^^ A(M2J 



(16) 



To the order shown, the mass-dependent chiral condensate still follows from eq. (15), after making 
the above substitution in the partition function (111), viz., 



lim 



(17) 



This is completely analogous to the case of dynamical quarks |2|, where instead one should replace 
the corresponding parameter fig by its modification due to (12). 



Let us now compare this result with what one gets in ordinary quenched chiral perturbation theory, 
which is valid in the opposite finite-volume regime where V ^ Then, to one-loop order 0, 110] 



5](m,; 



1 + 



1 



NrF'^ 



aA{Ml) - (ml - aM^dMiMM; 



+ ... 



(18) 



which also immediately follows from the quenched Feynman rules given above. As in the case of 
dynamical fermions Q there is an overlapping region where the two different expansions match on 
to each other. To see this, let us evaluate the leading-order contribution to the laige-fi^ expansion 
based on the relations (|l^ and (|l5| ) . We thus need the asymptotic expansion of the partition function 
(p!^) for mass-degenerate fermions. To the best of our knowledge, such an expansion has not been 
worked out yet. However, we can make use of the known asymptotic expansions in sectors of fixed 
topological gauge field charge and then perform the sum over topology subsequently pT[| . We 

first note that the inclusion of an arbitrary vacuum angle 6 in the effective partition function (11^ ) can 
be re-absorbed into a shifted Sq, so that 



dUo exp 



-M,TV (Uo + U'q 



2Nr 



9F/V2f 



(19) 



The projection on sectors of fixed topological charge v is then performed by means of a Fourier 
transform Q. The integral is Gaussian, and the result is, for large four- volumes V, 



1 

2^ 



2-K 



1 



dUo (det UqY exp 



.Tr (Uo 



U 



(20) 



up to an irrelevant z^-independent normalization factor. We have then simultaneously identified the 
topological susceptibility |jlO[ 

^ ^ . (21) 



2Nr 



We now insert the large-;U„ series expansion oi Zi^{fiy) (up to /Xj,-independent factors that are irrelevant 
right here) M, |, 



(22) 
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into the definition of the quenched chiral condensate 



lim — — — In Z ( /i,, I = 1 H + 

When next we wish to sum over topology 



1 + 



1/4 



2^2 



+ 



(23) 



(24) 



we are faced with the fact that i^v^) ~ 1^, while ^t, here is kept finite as y — > oo. This means that 
it is not immediately obvious if there is a region in which we can expect to find overlap between the 
two perturbative expansions discussed abovej^ The precise requirement for the z^-average of the above 
first-order term in the expansion to be small is 



V > 



2NrmlY? ' 



(25) 



while we also need 



V « 



4m2S2 



(26) 



oo 



This can be achieved by having itlq/Nc <C F"^ /2, which formally can be met by sending 
while keeping all other quantities fixed^. As discussed above, this limit is within the region of validity 
of the finite-volume e-expansion, and it gives a consistent regime in which we should have simulta- 
neous validity of both types of expansions. We can then insert the relation (|2l|) into the I'-averaged 
condensate 



lim (^cr-^ 



1 + ^ — - + 



1 + 



NrF'^M^V 



+ 



(27) 



and hence, from eqs. (|l7|) and (16), to leading order 



1 + 
1 + 



rrin 



NcF'^M'^V 

^ «A(M2 



1 + 



1 



aA(M, 



vvJ 



ml 



NrF"^ 



NcF^ 

{ml - aMDdMiMM, 



aMl)dM2^A{Ml) 
+ ... . 



+ 



(28) 



Thus, in this particular regime the leading-order correction to the chiral condensate precisely coincides 
with that of one-loop chiral perturbation theory (IH 



An alternative way of interpreting the formulas (16) and ( p^ is that, effectively, and to leading order 
in the finite- volume perturbation theory, one is working with a shifted quenched chiral condensate at 
finite volume: 



'eff 



1 



NcF^ 



aA{Ml) - (ml - aM^dMiMM, 



(29) 



^Indeed, the leading-order term in the expansion, which ordinarily is responsible for an overlapping region of validity 
1^, ^ of the two expansions is proportional to A^„, and hence vanishes in the quenched case. As shown below, in the 
quenched theory the region of agreement is just pushed one order higher. 

•^In a proper scaling with Nc, mo is of order 1, while is of order Nc, but this only makes it more clear that the 
inequality can be met. 
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It is clearly of interest to have this volume dependence made more explicit. So far the calculation has 
been made at fixed ultraviolet cut-off A[/y (having in mind a lattice cut-off with scale Af/y ~ 1/a). The 
explicit volume dependence is recovered once one notices that in most regularization schemes the boson 
propagator at the origin (^0|), and any derivative thereof, separates into ultraviolet divergent volume- 
independent terms, and additional finite volume-dependent terms [12|. The theory is renormalized as 
at infinite volume, and the explicit volume dependence is what remains. A detailed analysis has been 
made by Hasenfratz and Leutwyler |12|, from where we can borrow the main results. In dimensional 
regularization, and setting the subtraction point equal to unity, one has 



167r2 



[ln(M2) + ci]+5i(M2,Li) 



(30) 



where ci is a mass and Li independent constant that diverges near four dimensions. We have ignored 
the tree level contribution from the O(p^) terms in the chiral Lagrangian; they give volume-independent 
terms that are of no interest here. All volume dependence is contained in the functions (in four 
dimensions; for A{M^) we need only r = 1) ||2|, |l^ 



1(1?' 
V \ 4tt 



/ dt exp 
Jo 



M^LH 



l[Sm/{tL))-l] , 



where 



S{x) = ^ exp[— vrn^x] 



(31) 



(32) 



The general relation 



9r{M\Li) 



gr-i{M\Li) 



(33) 



allows one to calculate all ^^'s from go{M^ , Li). For small masses a systematic expansion has been 



given in ref. [12 



go{M^,Li 



|mml) + ^m^ 



In(ML) - - 



1 °° 1 
+ TTT.-ML^/L)iMLf- 



n=0 



(34) 



where (3n{Li/L) are "shape coefficients" that depend on the shape of the four- volume, and which can 
easily be evaluated numerically using the formulas of ref. |l^. Subtracting the zero-momentum pole 
term, we thus have 



9i{M\Li) 



gi{M\Li) 
1 

8^ 



1 



In(ML) - 



n=l 



'n-l)\ 



f3n{Li/L)M'^^'^-'^'^L^''-^ 



and 

UM\L,) = 
from which we obtain both 
A(m2) = 



1 -, oo 

ln( ML) + ^ + y 

n=2 



(n-2) 



■pn{Li/L)M 



2{n-2)^2n-4 



/l/f2 °° 1 

n=l ^ ' 



l3n{Li/L)M'^^''-^'^L^''-^ 



(35) 
(36) 

(37) 
(38) 
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and its derivative in a small-mass expansion. All volume dependence of the one-loop result (|29| ) is 
thus explicitly given. 



We note two important differences between the one-loop correction to the quenched chiral condensate, 
and that of the full theory with dynamical quarks: 



• The expression for the effective finite-volume chiral condensate ( |29D involves two new constants 
(mo and a), in addition to the pion decay constant F. 



• In contrast to the case of dynamical quarks (|12D, the leading correction contains not just the 
finite volume function gi{M^^, Li), but also^2(-^^^, -^i)- 

The first point means that it is not straightforward to use the quenched chiral condensate to extract 
the (quenched) pion decay constant F. The second point implies, as we will see shortly, that the 
finite-volume perturbative expansion of this quantity contains a logarithmic divergence. 

We wish to find the effective chiral condensate at the finite volumes considered here. To avoid spec- 
ifying the scale entering the logarithm and all volume- independent terms, it is convenient to simply 
compare Se// at two different four- volumes Vi and V2 of same geometry (here, let Li^2 = ^1^2^) • Based 



on (p9) and (pq), we find 



{V2) N,F^ 



+ ... (39) 



after having taken the limit My„ 0. In contrast to the case of dynamical quarks, what should be 
a small finite-volume correction can in fact become arbitrarily large. This logarithmic divergence of 
the effective chiral condensate as the volume V is sent to infinity is a direct analog of the well-known 
quenched chiral logarithms ||5|, ^ of infinite-volume quenched chiral perturbation theory. We call it a 
quenched finite volume logarithm. 

A logarithmic divergence in what should be a small correction of course signals that this finite-volume 
perturbative expansion of the quenched effective theory is problematic. The trouble that shows up as 
ln(Mt,t,) for ^ in the regime L ^ l/M^^, shows up as ln(L) in the finite-volume regime where 
L <C Roughly speaking, as l/M^y is lowered and crosses L, the infrared cut-off switches from 

being given by l/M^y to being given by L. For both expansions it is the very starting point of the 
expansion that appears internally inconsistent. It seems that at best one can hope that either of the 
two perturbative expansions have windows of validity where the one-loop corrections remain small. If 
we adopt this point of view, then we should clearly demand in our case that when comparing ^eff at 
two different volumes Vi and V2 we must require 

2 

'ln(Li/L2)| « 1 (40) 



in order to apply (39). It has been argued by Sharpe ^ that in usual chiral perturbation theory it 
may be possible to resum a certain class of diagrams that give rise to quenched chiral logarithms. 
However, also this partial resummation remains divergent. Presumably a similar resummation may 
be performed for the quenched finite volume logarithms. But we prefer to avoid speculation beyond 
the one-loop order to which we have computed the correction, and insist at least on the condition (pO|). 



*The function g2{M^g, Li) appears at next order in the unquenched case, but then it is suppressed by an additional 
factor of the volume V, and hence harmless. 
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Difficulties with quenched chiral perturbation theory in finite volumes (but associated with unusual 
power-law behavior in L) have been noted earlier in connection with the quenched pion scattering 
length 0]. 

A divergent chiral condensate for the quenched theory has long been suspected from many different 
arguments, of which the presence of ordinary (infinite-volume) quenched chiral logarithms ^ is the 
one closest to the present perspective. The divergence has been difficult to see in lattice gauge theory 
simulations, but a recent careful Monte Carlo study of the quenched Schwinger model using overlap 



fermions has revealed a clear increase in S as the volume V is made larger [14]. This should be directly 



related to the quenched finite volume logarithm found here, although the analysis will differ in detail 



due to the particularities of two space-time dimensions [15, 16] 



3 The Chiral Condensate and F from Dirac Operator Eigenvalues 



The dominance of the zero momentum mode in the effective chiral Lagrangian when l/M^^ ^ V is 
known to have a very direct consequence in terms the Dirac operator spectrum of the underlying theory 
R 17, |l8|. In this finite-volume range the spectral properties of the Dirac operator are analytically 



computable for an infinite sequence of eigenvalues. While originally derived on the basis of universal 



results from Random Matrix Theory [17, 18[, the low-energy spectrum of the Dirac operator can 
equally well be obtained directly from the effective Lagrangian framework [IC, 19, 20, 21 1. Let us 



denote the spectral density of the Dirac operator by p{X;mi, . . . ,miq^), where A is the eigenvalue. It 
was observed some time ago p3[ that the so-called microscopic spectral density WT 



1 (_^.J^ 



(41) 



gives an excellent finite-size scaling function from which to extract the infinite- volume chiral condensate 
S from finite lattice volumes. Even more information can be obtained from the individual eigenvalue 
distributions, which all depend on just this one single parameter S, and whose complete analytical 
expressions are now known [^] . Each single eigenvalue distribution can then be used to measure the 
infinite- volume chiral condensate S. 

The preceding section clearly suggests that apart from the chiral condensate one more quantity can be 
extracted from the microscopic Dirac operator spectrum, namely the pion decay constant F. While 
the quenched case is special because of the quenched finite volume logarithm, the situation is much 
more favorable when considering a theory with dynamical quarks. We recall that the effective partition 
function, including its one- loop correction (0), can be written 



-eff 



SU(Nf) 



dU exp 



(42) 



where 



1 



1 1 



Nf F2 

Taking the massless limit of eq. (^3|), the only correction to S comes from the finite four- volume: 

1 1 Pi{U/L) 



(43) 



'eff 



iV2 



1 + 



f 

Nf F2 L2 



+ 0(1/L^ 



(44) 
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For symmetric four-tori f3i = 0.140461... |12|, and for other shapes the coefficient (3i is readily computed 
numerically from the general formula of that same reference. The correction to S is thus explicitly 
known to one-loop accuracy, and it depends only on F. 



Prom the effective partition function (p2D, and its supersymmetric extension 19, 2C], follows 
the microscopic Dirac operator spectrum. To first order in perturbation theory it is unchanged in 
form compared with the leading behavior, except for the finite- volume shift indicated in (|4^). It is 
immediately clear from (^) that all predictions for fixed sectors of topological charge u are unchanged 
as well, except for the shift in S. With high enough statistics, the finite volume correction to T, can 
be measured from fits to different lattice volumes, yielding F. The two-loop contribution to the chiral 



condensate in this volume range has also been computed |24|, but at that order one can no longer 
view the correction as just a shift in S in an otherwise unaffected effective theory. The suggestion 
that the pion decay constant can be extracted using finite- volume scaling in this form was first made 



by Gasser and Leutwyler P], and in a slightly different setting by Neuberger [25|. The Dirac operator 



eigenvalues just give a particularly convenient way to make use of the formalism. 

The possibility of measuring the pion decay constant F by considering the combination of chiral 
perturbation theory and the analytical results for the microscopic Dirac operator spectrum was first 
explored by Berbenni-Bitsch et al. in ref . . They compared lattice Monte Carlo data for the chiral 
condensate in quenched SU(2) gauge theory with a general expansion of for staggered fermions 

at finite four-volumes. Relating one of the coefficients in the expansion to the pion decay constant F 
by means of the Gell-Mann-Oakes-Renner relation, they were thus able to extract F from a general 
fit of S(m) over a mass range of several orders of magnitude. A similar analysis has very recently 
been made for staggered fermions and gauge group SU(3) |2^. In that case there are independent 
estimates of F from earlier lattice studies, and the authors of ref. find very good agreement. 



A comment on staggered fermions can be made here. Comparisons between the analytical prediction 
for the Dirac operator spectrum of dynamical staggered fermions will typically be made so far from 
the continuum that the continuum number of flavors is not yet seen. Instead, the theory appears for 
just one staggered fermion to be of Nf = 1 |28], and in the effective theory the remaining bosons are 
heavy on the scale of the one single pseudo-Goldstone boson. The symmetry breaking pattern is thus 
U(l) xU(l)^U(l), and it is of interest to know also the leading correction for an unquenched theory 
with this Abelian Goldstone manifold. Doing the one-loop flnite-volume perturbative calculation we 
immediately find the analog of eq. (|l2|) : 

/ '^^ ^^)^ = + U^o)^iMl) , (45) 



where in this case Uq £ U(l). The effective chiral condensate is then rescaled according to 



^eff 



1 - ^^(0) 



(46) 



Thus also in this Nf = 1 theory one can extract the pion decay constant F from the smallest Dirac 
operator eigenvalue distributions. The general case of Nf staggered fermions at strong coupling follows 
similarly (a factor of A'^^ in front of the correction term in (4(:)). 



The above result for the Nf = 1 theory with staggered fermions away from the continuum limit may 
appear to contradict a recent lattice Monte Carlo study of the distributions of the smallest Dirac 
operator eigenvalues in SU(3) lattice gauge theory with one staggered fermion ||2^. There a very 
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accurate agreement was found when comparing with the analytical prediction based on the continuum 
Nf = 1 theory, where there is no Goldstone manifold at all, and all mesons are massive even in the 
chiral limit. In fact, there is no contradiction. The agreement found in was with the sector of just 
topological gauge field charge u = 0. The projection on this sector is actually equivalent to having 
a U(l) Goldstone manifold, as needed in order to compare with the Nf = 1 theory with staggered 
fermions away from the continuum. In detail, the projection on the z/ = sector, 

Z^=o = — r^de Z{e) = I dU e'^™"^^) (47) 
2vr 7o Ju{i) 

is completely equivalent to the zero momentum mode effective partition function for the coset of 
U(l)xU(l)^U(l) spontaneous symmetry breaking, without any projection on the z/ = sector. The 
role of the zero-mode U(l) field is simply being played by the angular integration that projects down 
on = 0. This also gives a complementary viewpoint to the observation that staggered fermions seem 
insensitive to gauge field topology at these couplings ||2^ (see ref. [^] for the picture that emerges as 
the continuum limit is taken). 



4 Conclusions 

We have analyzed quenched chiral perturbation theory in the finite-volume regime where the pseudo- 
Goldstone masses have Compton wavelengths much larger than the linear extent of the volume, 
» V. In this regime a good starting point for perturbation theory is not the usual plane- 
wave solution of the fields. Rather, it is a collective variable which is chosen to be the zero momentum 
mode of the given field. Using this perturbation theory, and adding the (uniterated) vertex corre- 
sponding to the singlet field ^o{x) as usual for the quenched theory, we have computed the one- loop 
correction to the quenched chiral condensate. It turns out that this correction becomes large when 
the volume is taken to infinity, in the chiral limit as a logarithm in the volume V. This so-called 
quenched finite volume logarithm prohibits taking the volume V to infinity. Whether or not it sig- 
nals a true divergence in the quenched chiral condensate clearly cannot be ascertained within this 
perturbative framework alone. Quenched finite volume logarithms are likely to be a generic feature of 
quenched chiral perturbation theory in the present finite-volume regime, and should thus appear in 
other observables as well. 

In the full theory with dynamical quarks the one-loop correction to the chiral condensate in this 
finite- volume regime amounts to a rescaling of S in an otherwise unchanged effective theory. We have 
used this fact to point out that the microscopic spectrum of the Dirac operator, and in particular the 
individual smallest eigenvalue distributions, can be used to extract the pion decay constant F. This 
entails determining numerically the leading correction to the well-known microscopic scaling for the 
Dirac operator spectrum near the origin, and it appears to be feasible. It is quite different from the 
program to extract these quantities and further couplings in the chiral Lagrangian from (ratios of) 
matrix elements whose expansions in partially quenched chiral perturbation theory are known |3l|, |3^ . 
So far the eigenvalue approach has only been tested away from the continuum and only for the case 
of S. It would be interesting to see if with present lattice volumes the finite- volume corrections can 
be determined to the accuracy required to extract F, and can be pursued towards the continuum. 
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